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Abstract 

In this work we study the global approximate multiplicative controllabil- 
ity for the linear degenerate parabolic Cauchy-Neumann problem 

v t - (a(x)v x ) x = a(t,x)v in Q T = (0,T) x (-1,1) 

a(x)v x (t,x)\ x=± i =0 t£(0,T) 

w(0, x) = vo(x) x e (— 1, 1) , 

with the bilinear control a(t,x) £ L°°(Qt)- The problem is strongly degen- 
^Nj erate in the sense that a G C ([—1,1]), positive on (—1,1), is allowed to 

CO vanish at ±1 provided that a certain integrability condition is fulfilled. We 

s ! will show that the above system can be steered in L 2 (Q.) from any nonzero, 

nonnegative initial state into any neighborhood of any desirable nonnegative 
\Q target-state by bilinear static controls. Moreover, we extend the above result 

f^) relaxing the sign constraint on vq. 
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1 Introduction 



Motivation 

Climate depends on various parameters such as temperature, humidity, wind in- 
tensity, the effect of greenhouse gases, and so on. It is also affected by a complex 
set of interactions in the atmosphere, oceans and continents, that involve physical, 
chemical, geological and biological processes. 
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One of the first attempts to model the effects of interaction between large ice 
masses and solar radiation on climate is the one due, independently, by Budyko 
0E] and Sellers [25] (see also [T2], HSl HZ] and the references therein). Such a model 
studies how extensive the climate response is to an event such as a sharp increase 
in greenhouse gases; in this case we talk about climate sensitivity. A process that 
changes climate sensitivity is called feedback. If the process increases the intensity 
of response we say that it has positive feedback, whereas it has negative feedback if 
it reduces the intensity of response. 

The Budyko- Sellers model studies the role played by continental and oceanic 
areas of ice on climate change. In such a model, the sea level mean zonally averaged 
temperature u(t, x) on the Earth, where t denotes time and x the sine of latitude, 
satisfies the following degenerate Cauchy-Neumann problem (1.1 1 in the bounded 
domain ( — 1, 1). 

The effect of solar radiation on climate can be summarized in the following 
figure 




Figure 1: www.edu-design-principles.org (copyrighted by DPD) 



We have the following energy balance : 



Heat variation = R a —R e +D 

• R a = absorbed energy 

• R e = emitted energy 

• D = diffusion 



The general formulation of the Budyko-Scllcrs model on a compact surface M. 
without boundary is as follows 



u t - A M u = R a (t, x, u) - R e (t, x, u) 
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where u(t, x) is the distribution of temperature and is the classical Laplace- 
Beltrami operator. Moreover, 

• R a (t,x,u) = Q(t,x)/3(x,u) 

• R e (t, x, u) = A(t, x) + B(t, x)u 

In the above, Q is the insolation function, and (3 is the coalbedo function (that is, 
I- albedo function). 

Albedo is the reflecting power of a surface. It is defined as the ratio of reflected 
radiation from the surface to incident radiation upon it. It may also be expressed 
as a percentage, and is measured on a scale from zero for no reflecting power of a 
perfectly black surface, to 1 for perfect reflection of a white surface. 




"Albedo" 

Ice/snow reflects 

Water absorbs 



Figure 2: www.esr.org (copyrighted by ESR) 

The main difference between Budyko's model and the one by Sellers, is that 
in the former the coalbedo function is discontinuous, while in the latter it is a 
continuous function. In fact we have 
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• Budyko 




f (3 u< -10 


f3(u) = < 


[/3 ,/3i] u = -10 




I Pi u > -10 , 


• Sellers 




' fto u < it- 


(3(u) = 


line it_ < It < W + 






where u± = — 10 ± S, S > 0. 



On .M = £ 2 the Laplacc-Beltrami operator is 

sin I 90 V <9</> / sin d\ 2 J 
where is the colatitude and A is the longitude. 




-50 -40 -30 -20 -10 10 20 30 



Annual Mean Temperature 

Figure 3: www.globalwarmingart.com (copyrighted by Global Warming Art) 



In the one-dimensional Budyko- Sellers we take the average of the temperature 
at x = cos (f> and the Budyko-Sellers model reduces to 

u t - ((1 - x 2 )u x ) x = g(t, x) h(x, u) + f(t, x), x E (-1, 1) 

(1.1) 

(1 - x 2 )u xlx=±1 = . 
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Problem formulation 

Let us consider the following Cauchy- Neumann strongly degenerate boundary lin- 
ear problem in divergence form, governed in the bounded domain (— 1, 1) by means 
of the bilinear control a(t, x) 

v t — {a(x)v x ) x = a(t, x)v in Q T = (0, T) x (-1, 1) 

a{x)v x {t,x)\ x=±1 =0 *G(0,T) (1.2) 

v(0, x) = vo(x) X £ (— 1, 1) . 

We assume that 

1. v eL 2 (-l,l) 

2. a£L°°(Q T ) 

3. a G C^fl-l.l]) satisfies 

(a) a{x) > Va: G (-1,1), a(-l) = a(l) = 

(b) A G L x (-1, 1), where A(z) = /* ^ . 



Remark We observe that 

1. A L x (— 1, 1), so a(x) is strongly degenerate 



2. the principal part of the operator in ( 1.2 ) coincides with that of the Budyko- 



Sellers model for a(x) — 1 — x 2 . In this case A(x) = |ln ( jz^j £ ^(—1, 1) 

3. a sufficient condition for 3.b) is that a'(±l) 7^ (if a G C 2 ([— 1, 1]) the above 
condition is also necessary). 

□ 



We are interested in studying the multiplicative controllability of problem ( 1.2 1 
by the bilinear control a(t,x). In particular, for the above linear problem, we will 
discuss results guaranteeing global nonnegative approximate controllability in large 
time (for multiplicative controllability see [5]). 
Now we recall one definition from control theory. 

Definition 1.2 



We say that the system (1.2) is nonnegatively globally approximately control- 
lable in L 2 {— 1,1), if for every e > and for every nonnegative vq(x), Vd(x) G 
L 2 {— 1, 1) with vq ^ there are a T = T(e,Vo,Vd) and a bilinear control a{t,x) G 



L°°(Qt) such that for the corresponding solution v(t,x) of (1.2) we obtain 

\\v(T,-)-v d \\ L 2 ( _ 1A) <e. 
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In the following, we will sometimes use || • || instead of || • 1). 
Main results 

In this work at first the nonnegative global approximate controllability result is 



obtained for the linear system (1.2 1 in the following theorem. 
Theorem 1.3 



The linear system (1.2) is nonnegatively approximately controllable in L - (— 1, 1) 



by means of static controls in L°°{—\, 1). Moreover, the corresponding solution to 



1.2) remains nonnegative at all times. 



Then the results present in Theorem |1.3| can be extended to a larger class of 
initial states. 

Theorem 1.4 

For any Vd € L 2 (—l, l),v<j > and any Vq € L 2 (— 1, 1) such that 

l 

VQVddx > 0, (1-3) 

l 

for every e > 0, there are T — T(s,VQ,Vd) > and a static bilinear control, 
a = ct(x), a € L°°{— 1, 1) such that 

\\v(T,-)-v d \\ L 2 ( _ 1A) <e. 



Remark The solution v(t, x) of the problem ( 1.2 ) in the assumptions of Theorem 
1.4| does not remain nonnegative in Qt, like in Theorem |1. 3 1 but it can also assume 
negative values. q 



Mathematical motivation 

This note is inspired by [501 [5]. I n EH] A.Y. Khapalov studied the global nonnega- 
tive approximate controllability of the one dimensional non- degenerate semilinear 
convection-diffusion-reaction equation governed in a bounded domain via the bi- 
linear control a € L°°(Qt)- In [8], the same approximate controllability property 
is derived in suitable classes of functions that change sign. 

In this note we extend some of the results of J2Uj to degenerate linear equations. 
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General references for multiplicative controllability are, e.g., [T5I UM I2T1 [251 15] . 
In control theory, boundary and interior locally distributed controls are usually 
employed (see, e.g., [H [101 HU HH HS1 US]). These controls are additive terms 
in the equation and have localized support. However, such models are unfit to 
study several interesting applied problems such as chemical reactions controlled 
by catalysts, and also smart materials, which are able to change their principal 
parameters under certain conditions. This explains the growing interest in multi- 
plicative controllability. 

2 Preliminaries 

Positive and negative part 

Given f2 C K n , v : f2 — > K we consider the positive-part function 

v + (x) = max (v(x), 0) , Vx G f2 , 

and the negative-part function 

ir(.i) =max(0,-'o(s)), Vx € £1 . 

Then we have the following equality 

v = v + — v~ in 

For the functions v + and v~ the following result of regularity in Sobolev's spaces 
will be useful (see [21], Appendix A ). 

Theorem 2.1 

Let n C R n , u-.fl — > R, u e H 1 ' 3 ^), 1 < s < oo. TTien 

and /or 1 < i < n 

!u x . in {x £ ft : u{x) > 0} 

(2.4) 
in {x E CI : u{x) < 0} , 

and 

-u Xi in {x £ ft : u{x) < 0} 

(«")*, = 4 (2-5) 
in {a: € Q : m(x) > 0} . 
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Gronwall's Lemma 



Lemma 2.2 Gronwall's inequality (differential form). 

Let r](t) be a nonnegative, absolutely continuous function on [0,T], which sat- 
isfies for a.e. t G [0,T] the differential inequality 

T)'(t)<<j>(t)r](t)+TP(t), (2.6) 



where <j)(t) and tp(t) are nonnegative, summable functions on [0,T]. 
Then 



n(t) < eti ^ ds 



ij(0)+ / V(s)ds 



(2.7) 



/or allO <t<T. 



In particular, if tp(t) = m (2.6), i.e. r)' < (j>T) for a.e. t £ [0, T], and 77(0) = 0, 

t? = m[0,T]. 



Well-posedness in weighted Sobolev spaces 



In order to deal with the well-posedness of problem (1.2 1, it is necessary to intro- 
duce the following Sobolev weighted spaces 

ffa(-M) : = 

:= {u e L 2 ( — 1, 1) : it locally absolutely continuous in (—1, 1), y/au x € L 2 (— 1, 1)} 
and 

ff a 2 (-l, 1) := G i$(-l, 1)| a Ul G H\-l, 1)} = 
= {ii£ L 2 ( — 1, locally absolutely continuous in (—1, 1), 

ra£ £$(-1,1), (Hi, G H^-ljl) and (au x )(±l) = 0} 
respectively with the following norms 

\\ u \\m a ■= \\ u \\h(-i,i) + ML and \\u\\ 2 H 2 := \\uf H t + \\(au x ) x f L2{ _ 1A y, 

where |tt|i )0 = ||v^ u :e||l 2 (-i,i) is a seminorm. 

In this note we obtain the following result. 

Lemma 2.3 

H a (— 1, 1) <—> L 2 (— 1,1) uiii/i compact embedding. (2-8) 
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Proof: 

Given u £ 1), let 



u(x) 



u if x e [— 1, 1] 

elsewere . 

It is sufficient to prove that, for every R > 0, 



sup / \u(x + h) — u{x)\ dx — > 0, 
|u||i,o<R. 



as /i -> (2.9) 



Let h > 0(Q and let u G 1, 1) be such that ||w||i, a < -R, we have the following 

equality 

|tt(x + /i) — u(x)\ 2 dx = 



-1 r-l-h 

\u(x + h)\ 2 dx + I \u(x + h) - u(x)\' 2 dx + / \u(x)\ 2 dx 

l-h 



l-h 
-1+h 

-1 



-1 

l-h 



\u{x)\ 2 dx + / \u(x + h) - u(x)\ 2 dx + / \u(x)\ 2 dx 



l-h 



First, let us prove that 



sup 

l|w||l,a<H»'-l 



\u(x + h) - u(x)\ 2 dx — >0, as/i^0 + . (2.10) 



Recalling that A(x) — J Q ^y, we have 



rx+h | 

\u(x + h) - u(x)\ < / y/g(s)\u' (s)\ ds < 

Jx V"( S ) 



\ 2 / l-X+h J 

// \|2 • > ' ' 



< / a{s)\u'(s)\ 2 ds 



a(s) 



\u\x.a [A(x + h) -A{x)Y 



By integrating on [—1, 1 — h], since A € L (— 1, 1) (by assumption 3.b)), we obtain 



i-/t 



|«(a; + h)- u(x)\ 2 dx < \u\ 2 l a / (A(x + h) - A{x)) dx < 



l-h 



< ir 



l+h 



l-h 



A(x) dx - I A(x) dx 



1 In the case h < we proceed similarly. 
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R 2 



A(x) dx 



1-h 



-1+h 



A(x) dx 



0, as h -> H 



Now, let us prove that 



sup / \u(x)\ 2 dx — > 0, 

\\u\\i, a <RJl-h 



as 



h-> 0+ 



(2.11) 



We have 



|«(0)| < |«(»)| + [ X yMs)\u'(s)\^t=ds< 
Jo v a ( s ) 



< \u(x)\ + (J i a( S )\u'(s)\ 2 ds^j 2 Qf \ 2 < !„(»)! + |„| 1>aV G4(i). 
By integrating on [0, 1], we obtain 

\u(0)\< [ \u(x)\dx+ \u\i, a [ \/A{x) dx < 



< ||u||z = (-l,l) + |u|i, / y/A(x)dx < C||u||i, a 
Jo 



Then, 

Now, it follows that 



|w(0)| < CR. 



(2.12) 



\u(x)\ 2 < 2\u(0)\ 2 + 2A(x)\u\l a <CR 2 + 2A{x)R 2 . 
Finally, since A 6 L (— 1, 1), by integrating on [1 — h, 1] we obtain 



/ \u(x)\ 2 dx < C hR 2 + 2R 2 [ A(x) dx — ► 0, 

Jl-h Jl-h 



as h -> H 



Similarly, we can prove that 
-l+ft 



sup / \ u i x )\ dx — > 0, 

\H\l, a <RJ-l 



as ft, -> 0^ 



(2.13) 



By (2.101, (2.11) and (2.13) we obtain (2.9) 



□ 



We now recall the existence and uniqueness result for system (1.2) obtained in 



[7] (see also [I]). Let us consider, first, the operator (A , D(A )) defined by 

D(A ) = H 2 (-l,l) 
A u = {au x ) x , VuG D(A ) . 



(2.14) 
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Observe that Aq is a closed, self-adjoint, dissipative operator with dense domain 
in L 2 (— 1, 1). Therefore, Aq is the infinitesimal generator of a Co — semigroup of 
contractions in L 2 (— 1, 1). 

Next, given a £ L°°(— 1, 1), let us introduce the operator 



(2.15) 



D(A) = D(Aq) 

A = A + al . 
For such an operator we have the following proposition. 

Proposition 2.4 



D(A) is compactly embedded and dense in L 2 (— 1, 1). 

A : D{A) — > L 2 (—l, 1) is the infinitesimal generator of a strongly continu- 
ous semigroup, e tA , of bounded linear operators on L 2 {— 1, 1). 



Observe that problem (1.2) can be recast in the Hilbert space L 2 (— 1, 1) as 

u'(t) = Au(t) , t > 

(2.16) 

u(0) = Uq . 
where A is the operator in ( |2.15[ ). 



We recall that a weak solution of (2.161 is a function u G C°([0, T\. L 2 (— 1, 1)) 



such that for every v £ D(A*) the function (u(t),v) is absolutely continuous on 
[0, T] and 

±{u(t),v) = {u(t),A*v), 
for almost t € [0,T] (see 0). 



Theorem 2.5 

For every a <E L°°((0,T) x (—1,1)) and every uq € L 2 (— 1,1), there exists 
unique 

u e C°([0,T}; L 2 (-l, 1)) n L 2 (0,T; ff^-l, 1)) 
weafc solution to (1-2), which coincides with e tA u$. 



In the space 

B(0, T) = C°([0, T];L 2 (-1, 1)) n L 2 (0,T; ff^-1, 1)) 
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let us define the following norm 

l-T f l 



||«Hil ( o ) D= sup ||u(t,-)||l a( _i,i)+2 / / a(x)uidx, Vu€B(0,T). (2.17) 
te[o,T] Jo J-i 

3 Some auxiliary lemmas and the proofs of main 
results 



Let A = Aq + al, where the operator Aq is defined in (2.14) and a <E L°°(— 1, 1) . 
Since A is self-adjoint and D(A) L 2 (— 1,1) is compact (see Proposition 2.4), 
we have the following (see also [4]). 

Lemma 3.1 

There exists an increasing sequence {Xk}keNi with — > +oo as k — > oo , 
such that the eigenvalues of A are given by {— Xk}keN, oxid the corresponding 
eigenf unctions {u)k}k&$ form a complete orthonormal system in L 2 (— 1, 1). 

In this note we obtain the following result 
Lemma 3.2 

Let v e <7°°([-l,l]),t; > on [-1,1], let o,(x) = x € (-1,1). 

Let A be the operator defined in 1 2. Iffy with a = a* 



D{A)=Hl(-lA) 

(3.18) 

A = A + ctj , 

and let {\k}, {&k} be the eigenvalues and eigenf unctions of A, respectively, given 

by Lemma\3.1\ Then 

v 

Ai = and \u>i\ = r. — rr. 

IMI 

Moreover, w^n and — ^ are the only eigenfunctions of A with norm 1 that do not 
change sign in (—1, 1). 



Remark Problem (3.18) is equivalent to the following differential problem 



(a(x)uj x ) x + a*(x)uj + \uj = Q in (—1,1) 

(3.19) 

a(x)ui x (x)\ x=±1 = . 

□ 
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Proof: (of Lemma 3.2 ) 
STEP.l Wc denote by 



and 



{^k}ki 



respectively, the eigenvalues and orthonormal eigenfunctions of the operator (3.18) 
(see Lemma 3.1). Therefore, 

(wfc,Wfc) L 2 ( _ M ) = J w k (x)w h {x)dx = 0, lih^k. (3.20) 

We can see, by easy calculations, that an eigenfunction of the operator defined in 
(3.18) is the function 

v(x) 



associated with the eigenvalue A = 0. Taking into account the above and consid- 
ering that v(x) > 0, Va; e (-1, 1) 

Ik* GN : Wfc » = ^>0 or cj kt (x) = - V ^4<0, Vie (-1,1). (3.21) 
M \\v\\ 



Writing (3.20) with k = fc* we obtain 



(wk tl W)i)iii(-i,i) = / u k *{x)uj h {x)dx = 0, V/i^fc» 



(3.22) 



Therefore, considering (3.22) and keeping in mind that uj^, > or ui k , < in 
(—1, 1), we observe that u) kt is the only eigenfunction of the operator defined in 
( 3.18[ ) that doesn't change sign in (—1,1). 



STEP.2 Let us now prove that 



k jf — 1 . 



(3.23) 



that is, Ai = 0. By a well-known variational characterization of the first eigenvalue, 
we have 

/, (flUj — a* u 2 ) dx 

Ai = inf ; . 

ueffi(-i.i) j\u 2 dx 



By Lemma 3.1 since A/^ = 0, it is sufficient to prove that Xi > 0, or 



a* u 2 dx < 



a u x dx, 



V w e^(-i,l) 



(3.24) 



Integrating by parts, we have 



a* u 2 dx = — 



1 {avx)x 
-l v 



u 2 dx = 



dx 
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2uu„ 



a v x dx 

-i v 



= 2 / \/a —u\fau x dx — ai 
J-i v J_ x 



< 



v x u 

V 



dx 



, dx 



dx < 



dx = 



dx . 



from which (3.24). 



□ 



For the proof of Theorem 1.3 the following Lemma is necessary. 
Lemma 3.4 

Let T > 0, a G L°°(Q T ), let v Q £ L 2 (-l, 1), v (x) > a.e. x G (-1, 1) and 
let v G 6(0, T) be the solution to the linear system 



v t - {a{x)v x ) x = a(t,x)v 
a(ar)i; a! (t,x)| a!= ±i = 
v(0,x) = v (x) 



Q T = (0,T) x (-1,1) 

te(o,T) 
xe (-1,1) . 



The 



v(t,x) > 0, V(t,x) G Qq 



Proof: Let v G 6(0, T) be the solution to the system (1.2), and we consider the 
positive-part and the negative-part. It is sufficient to prove that 

v~ (t, x) = in Qt ■ 



Multiplying both members equation of the problem ( 1.2 1 by v and integrating it 
on (— 1, 1) we obtain 



[v t v — (a(x)v x ) x v — aw ] dx = 0. 



(3.25) 



Recalling the definition v + and v , we obtain 



VtV dx — I (v + — v ) t v dx = — / (v ) t v dx 



1 d 

~2~dt 



(v'fdx. 



Integrating by parts and applying Theorem 2.1 we obtain v G 1, 1) and 

the following equality 

/ (a(x)v x ) x v~ dx — [a(x)v x v^] 1 L 1 — / a(x)v x (— v) x dx = / a(x)v x dx. 
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We also have 



avv dx 



a(v ) 2 dx 



and therefore (3.251 becomes 

1 d " X 



from which 



~2Jt ' ( v ~> dx+ I a ( v ~) dx 



a{x)vl > 0, 



(v~) 2 dx < 2 / a(v~) 2 dx < 2||a|| 00 / (v~) 2 dx 



dt 

From the above inequality, applying Gronwall's lemma we obtain 
J (v~(t,x)) 2 dx<e 2tM ™ J {v-(0,x)) 2 dx. 

Since 
we have 
Therefore, 

v~(t,x) = 0, V(t,x)eQ T . 
From this, as we mentioned initially, it follows that 

v(t,x) = v + (t,x) > V(t,i)eQ T . 
We are now ready to prove our main result. 



v(0, x) = vq(x) > , 
v-(0,x) = 0. 



□ 



Proof: (of Theorem 1.3 ) 

STEP.l To prove Theorem |1.3| it is sufficient to consider the set of target states 
^eC°°([-l,l]), » d >0on[-l,l]. (3.26) 

Indeed, regularizing by convolution, every function Vd G L 2 {— 1, l),Vd > can be 
approximated by a sequence of strictly positive C°°([— 1, 1])— functions. 

STEP.2 Taking any nonzero, nonnegative initial state vq G L 2 {— 1,1) and any 
target state Vd as described in (3.26) in STEP.l, let us set 



a* (a;) 



(a(x)v dx (x)) x 
Vd(x) 



xe (-1,1). 



(3.27) 
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Then, by (3.26), 



We denote by 



{~A fc } 



a. (as) G L°°(-l,l) 



and 



A-ef 



respectively, the eigenvalues and orthonormal eigenfunction^of the spectral prob- 
lem Aui + \lj = 0, with A — Aq + a*I (see Lemma 3.1 ). 



We can see, by Lemma [3T2| that 
Ai = and 



Ul (x) = ^4 >0, Vz e (-1,1). 



(3.28) 



STEP.3 Let us now choose the following static bilinear control 

a(x) = a*(x) + (3, Vx G (—1, 1), with /3 G K (/? to be determined below). 



The corresponding solution of (1.2), for this particular bilinear coefficient a, has 
the following Fourier series representation Q 

v(t,x) = £y-A*+fl*^y ^ Q ( s )u, fe ( s )^o; fc (x) = 

= ( y ^ V*M(*)^i(z) + £ e <-**+«* ( y ^ Wo ( s H( s )d S ^ fc (x) 

Let . 

r(t,x) = X;e ( - A * +Wt ( / i*(a)w fc (a)dsWaO 
fe>i \J-i J 

where, recalling that X k < \k+i, we obtain 

-A fc < -Ai = for ever k G N, fc > 1 . 



Owing to (3.28), 



K*,-)-«*ll < 



v (s)uji(s)ds wi - 



|r(t,a;)||= 



\r(t,x)\\ 



2 As first eigenfunction we take the one which is positive in (—1, 1). 

3 Observe that adding /3 £ R in the coefficient c*» there is a shift of the eigenvalues corre- 
sponding to a t from {— AfcjfcgN to {— + /3]-fceN > but the eigenfunctions remain the same for 
a* and c?„ + /3. 
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Since — Afe < — A2, Vfc > 2, applying Parseval's equality we have 

2 



\\r(t,x)f<e 2 ^ + ^Yl f Wkda 



||wjt(x)|| 



= 2 (-A 2+ « t ^ (uo ^ fc)2 = e2 (-A 2+W||wo|| 2_ 
k>l 



Fixed e > 0, we choose T e > such that 



e -A 2 T e = £ 



f\ v v d dx 

WINP 



(3.29) 



Since wq € L ( — 1, 1), > and ^ in (—1, 1) and by (3.28), we obtain 



(«o,wi) 



Then, it is possible choose /3 £ so that 



i>o(a;)a;i(x)c?a; > 0. 



(3.30) 



VQUidx = \\v d \\ , 



that is, since wi = tt^-tt • 
' 1 IMI ■ 



IMI 2 \ 
Te \j\voVddxJ 



In 



(3.31) 



So, by (3.29), (3.31) and the above estimates for \\v(T £ ,-) - v d (-)\\ and ||r(T e ,-)|| 



we conclude that 



\\v(T e ,.)-v d (.)\\ <e<--*+W>\\M =e-^ ^ \\v \\ 

J_! v v d dx 



From which we have the conclusion. 



□ 



Proof: (of Theorem 1.4 1 
The proof of Theorem 
in STEP.3, inequality 1 



1.3 can be adapted to Theorem |1.4[ keeping in mind that 



3.30) continues to hold in this new setting. In fact we have 

v d (x) 



1 f 1 , 

v (x)oj!(x)dx = / v (x) 

1 J-i IMI 



dx 



1 



Ml J-i 



VQV d dx > 0, by assumptions (1.3) 



From this point on, one can proceed as in the proof of Theorem 1.3 
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